Representations in Density Dependent Hadronic Field Theory 
and compatibility with QCD sum-rules 

R. AguirrJ*! 

Departamento de Fisica, 
Facultad de Ciencias Exactas, 
Universidad Nacional de La Plata. 
C. C. 67 (1900) La Plata, 
Republica Argentina. 
(Dated: February 9, 2008) 

Abstract 

Different representations of an effective, covariant theory of the hadronic interaction are exam- 
ined. For this purpose we have introduced nucleon-meson vertices parametrized in terms of scalar 
combinations of hadronic fields, extending the conceptual frame of the Density Dependent Hadronic 
Field Theory. Nuclear matter properties at zero temperature are examined in the Mean Field Ap- 
proximation, including the equation of state, the Landau parameters, and collective modes. The 
treatment of isospin channels in terms of QCD sum rules inputs is outlined. 



* Electronic address: |aguirre@venus.fisica.unlp.edu.ar| 



1 



I. INTRODUCTION 



A generalized concept in present theoretical physics is that description of all physical 
phenomena should be derivable from first principles in a unified way. However, the bridge 
towards concrete applications requires elaborated procedures and judicious arguments. This 
is the case of Quantum Chromodynamics (QCD), which is the accepted theoretical model 
for the strong interactions. Despite the fact that it is perturbative in the high energy regime, 
the fundamental state of matter corresponds to the opposite limit, where confinement and 
the breakdown of symmetries make it mathematically intractable. Different effective models, 
such as Nambu- Jona Lasinio, Skyrme, bag-like, and chiral perturbation theory attempts 
to translate the main features of QCD into the hadronic phase. On the other hand, lattice 
simulations and QCD sum-rules share this aim, although using different methods. 
QCD sum-rules is a ingenious procedure to reveal the foundations of certain hadronic prop- 
erties. The method is based on the evaluation of correlation functions in terms of quarks and 
gluons degrees, then applying the operator product expansion it is possible to express it as 
combinations of perturbative contributions and condensates (non-perturbative). Finally this 
expansion is connected with the hadronic counterpart by means of a Borel transformation. 
The method was developed to study meson [l| as well as baryon j2]_properties in vacuum, 
it was subsequently generalized to study finite density systems js], \^. These calculations 
provide a useful guide for some static properties of hadrons immersed in a dense medium, 
as concerning QCD symmetries and phenomenology. However, they are not able to take 
into account the dynamical aspects of hadronic matter, which may be achieved by inserting 
coherently these results into a theoretical model of hadronic interactions. 
A similar situation was found in the study of nuclear structure, as the standard combination 
of microscopic potentials and the relativistic Dirac-Brueckner approximation for nuclear 
matter gives rise to reliable, density dependent, nucleon self-energies. Notwithstanding, this 
procedure was inadequate to treat finite nuclei due to unsurmountable mathematical diffi- 
culties. A feasible solution to this dilemma was proposed in 0,0, ^, by defining density 
dependent meson-nucleon vertices in terms of the self-energies obtained in Dirac-Brueckner 
calculations. The enlarged hadronic model keeps the mathematical versatility of the Quan- 
tum Hadrodynamic models but is equipped with couplings reflecting the properties of 
the nuclear environment. The structure of spherical nuclei, in the medium to heavy range. 
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was studied within this framework using the Hartree approximation. 

Further improvements, developed in , replace the density dependence of the meson-nucleon 
vertices by an expansion in terms of in-medium nucleon condensates. This conceptual re- 
placement restores the covariance and the thermodynamical consistency of the model, giving 
rise to the Density Dependent Hadronic Field Theory (DDHFT). It was recently appended 
by introducing momentum dependent vertices [l^, and a expansion of the vertices in terms 



of meson mean field values 



The validity of the method is only justified a posteriori, and relies in the flexibility of 
hadronic field models to accommodate pieces of information provided by another field of 
research. The enlarged model yields a simpler and more intuitive description, instead of 
making involved calculations based on first principles interactions. A significative exempli- 
fication of this standpoint was given by the Brown-Rho scaling of hadronic masses. Taking 
into account the chiral and scale symmetries of QCD solely, an approximate scaling law for 
the in-medium hadronic masses was derived in • This hypothesis was applied to describe 
heavy ion collision, reaching an excellent agreement with the experimental results for the 
low mass dilepton production rate Q|. 

Given a set of physically meaningful self-energies as input for DDHFT, there is room for 
a full family of hadronic models, according to the field parameterization assumed for the 
vertices. In principle, this may lead to different predictions for nuclear observables. This 
point has not been investigated yet, and it is the main purpose of the present work. 
In a previous paper , the author attempted to make use of the scheme outlined above to 
relate nuclear observables with QCD inspired results. For this purpose we took as input the 
nucleon self-energies in symmetric nuclear matter obtained in |4,] by using QCD sum rules. 



However, Ref. 



ll| does not exhaust the physical description given there, as it covers isospin 



asymmetric nuclear matter also. So, we aim to complete our theoretical model including 
isospin degrees of freedom, based on the results presented in 

We have organized this work presenting in the following section the general features of 
DDHFT, and particularly the meson parameterization of the hadronic vertices. In Section 
im we present and discuss the results for symmetric nuclear matter. A final summary is 
given in Section HVl 
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II. DENSITY DEPENDENT HADRONIC FIELD THEORY 



As mentioned above, DDHFT was designed to incorporate into the hadronic field formal- 
ism pieces of information produced by other theoretical frames. Strictly speaking, DDHFT 
takes as input the nucleon self-energies, which can be decomposed into iso-scalar and iso- 
vector components, each one containing Lorentz scalar and vector contributions: i. e. 
S^", S*"**, S^"", where the superscript a distinguishes iso-vector quantities, and a = 1,2 
stands for proton and neutron respectively. For homogeneous nuclear matter in steady 
state, they can be parameterized as functions of continuous parameters, like baryonic cur- 
rent density and temperature, characterizing the macroscopic state of hadronic matter. 
In particular, there exists a reference frame where the spatial part of the baryonic current 
vanishes in the mean. In its former version S*" was taken from relativistic Brueckner- 
Hartree-Fock calculations with one boson exchange potentials 

On the other hand, a lagrangian density is proposed in terms of meson fields a, uj^ in 
the iso-scalar sector, and (pc, p)^ the iso-vector one. The indices c = 1,2,3 stand for the 
meson isospin projection 

+ \{d,ad'^a - mW) - \w,.W^'' + \mluj,uj^ (1) 

where W^u = d^oj^—dyUfj,, R^^^^ = d^pl—d^p^^ has been used, and \I' represents a isospinor. We 
have adopted the convention that repeated isospin indices must be summed. The isovector 
vertices F*^ take values over the space generated by the Pauli matrices r and the identity. 
The masses have been fixed at the phenomenological values M = 940 MeV, = 783 
MeV, rris = 550 MeV, rrir = 770 MeV, and mj = 984 MeV. The iso-scalar vertices are 
assumed to depend on the scalar combinations: Si = ^JJ^-, Si = ■, ra^ = ^uj^uj^^ and 
m2 = 0", whereas the isovector F^ could depend linearly on jc = ^t'^'^ ^ or 0c, and on the 
scalars S3 = \J j'j^jc-, S4 = \fjc3c-, "^z = \j P^lPc-, and irii = \/(pc<Pc- The baryonic and isospin 
current have been written as = ^'7^\I' and = ^'7^r^\l/, respectively. This is not the 
most general dependence, but it keeps a clear separation between isospin scalar and vector 
degrees of freedom. In Ref. the dependence on si or S2 was presented and the case with 



only S2 was explicitly studied; in 



the meson dependence of the vertices was introduced 



within DDHFT, for symmetric nuclear matter. 

The field equations obtained for this lagrangian density are 



OlpA 



(□ + ml) a = ^ 
d^W^'' + mluj'' = ^ 
(□ + mj) 0" = 
d^R^" + my, = 



dm' 



-a 



dmn 



^ dmi ( dVs 



dm^ f dV) 
dmA 



dV^ 

dujy \dmi dmi 
dVt 



dmA 



dm, fdTJ dT^, 



(2) 
(3) 
(4) 
(5) 
(6) 



dp'l \dm?j^"' ' dm-, 

where the isospin indices have been denoted as: A, B = 1,2 for nucleons, and a, c = 1,2,3 
for mesons. In the first equation, the vertex derivatives must be interpreted as 
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dV dT dsk 

dip A dsk dip A ' 



The mean field approximation (MFA) is suited to describe homogeneous matter. Within 
this scheme, meson fields are replaced by their uniform mean values and bilinear combina- 
tions of spinors are replaced by their expectation values. For homogeneous, isotropic nuclear 
matter in steady state, some simplifications arise in the meson mean values. For instance, 
their coordinate dependence can be neglected; the spatial components of the vectors cu^ and 
Pfj, become null, and only the third component of the isovector mesons are non-zero, due to 
isospin conservation. 

Therefore within MFA, Eqs. (j21ISl) reduce to 



0= [(^^-M)5^B + S^^-^^^] iPb, 



,2- 



mlCu =< ^ T^Y"^ > - < ^' 



dm2 
^ dr, ^ auo' 



dmi mi 



^^^^ 
dm^ m^ 

dm, m, 
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(7) 


dm2 


(8) 


dT^ ^ u'^u'' 

-— 7^^^ >; 

dmi 


(9) 




(10) 


dm, m. 


(11) 



where tildes over the meson symbols stand for their mean field values, and the two last 
terms within brackets in Eq. ((7j) can be regarded as the scalar and vector components of 
the nucleon self-energies: 

Sf^ = < r,(T + ^A2^ >6aB+< rf + T^"" M4^^^ >, (12) 

S;^^ = < T^uj, - ^rAi^v]/ >6ab-< Tfj^pl + r^"" ^K^^^-^ >, (13) 

with Kk = d [Vsa - ^ + r^0, + ^,) /dsu . 

For the sake of concreteness, we have examined three possible parameterizations: 

a) symmetric nuclear matter, with rs(si), r^(si); 

b) symmetric nuclear matter, with rs(s2), r^(si); 

c) asymmetric nuclear matter, with Vs{jn2), r^(mi), rj(0'^,m4), = Fr(0'^,m3). 

The first two cases are comparable to the instances presented in [3] as vector and scalar 
dependencies. The last case is a development of the preliminaries calculations of 0, the 
field dependence of the isospin vertices have been chosen so as to adjust the QCD sum rules 
calculations of j^. 

These results are simplified within the assumptions (a-c) above, for the sake of completeness 
we consider separately each of these cases: 



a) = F^a, = g 



r^u; + - — ujub - - — aus 



(14) 



a = TsUs/m^, Co = T^nB/m^, 



b) = (y [Ts + '^nA ,Y.^ = g^QUJ [v^+'^^nB] , (15) 



aus J \ auB 

a, and tuhave the same expressions as in (a), 

„\ s^AB v-iisosr I \^isov AB s^AB „ /'visosr i ^isov AB\ t'tci\ 
C) = ^AB + ^ = 9u0 [^v ^AB + ^„ j , (16) 

Sisos "p ~ sr^isov AB -pSAi? J sr^isos -p ~ \^isov AB -pS AS ~ 

s — J- sCr, 2^^ — i J 9-, — w^-, — ~'- r Pi 

m^a = Us, m^uj = ub, 

da duo 

^ rf0 d<P dp 

Obtaining the equations listed above, we have used the Wick theorem and the Hartree 
approximation, therefore contractions of fields belonging to different full contracted terms 
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(both in Lorentz and isospin indices) have been neglected. This approach allowed us to 
extract the vertices F and its derivatives from the expectation values, and to express the 
mean value of products of more than two fermion fields as product of baryonic densities 
and currents. Finally, we have adopted the reference frame of static matter, therefore we 
have < jf, >= negf^o, < >= (n^J^ - n§)g^,Q6^^, with =< ipAlo'ipA > the number 
of neutrons or protons per unit volume, for A = 2,1 respectively, and ns = n^s^ + rS^ 
the baryonic number density. Furthermore we have introduced the simplifying notation 

ni^^ =< ipA'^^A >, ns =< ^'^ >, F =< F >, Si = Hb, S2 = n^, S3 =| n^^ - n^-* |, 

I (1) (2) I ~ ~ -A 1 

S4 =\ ris — Us I, mi = tu, m2 = 0", ms = p, and 1714 = (p. 

Cases (a) and (b) yield meson equations resembling those of the Walecka model Q], 

and nucleon self-energies with the rearrangement contributions added. In the instance (c) 

the self-energies have the same structure as in QHD-I, but the source term of the meson 

equations are modified. With respect to the isovector meson field equations, it must be 

noted that < ^Tc"^ >= for c = 1,2 as isospin is a conserved charge. In order to give 

explicit expressions for the corresponding vertices, we follow the parameterization of jj]: 

^AA{in) ak + IsPkt J . s 

= -^ -, 72 — 7) k=S,V (17) 

with 1^ = 1(— 1), for protons (neutrons), no = 0.17/m~^, t = [n^B^ —n^B^)/nB, the remaining 
coefficients can be expressed in terms of the quantities used in Ref. ^1] as A = C^VN + C^gN + 
CI, + a = Cl2 + CIA% a, = MX- Ci^^ - ClA^ = Ma - {^Cr, + CiAi]), 

a, = -{CPvN + MC^^ + MCPAfg), (3, = -{C^_v]^ + MC^^ + MCpAII), which were obtained 
in ^, by averaging over the Borel mass. For the sake of completeness we give here the 
numerical values vn = 3, f ^ = 1, gN = — 8M/9, kn = 8, (n = 0.54, Cj. = —0.042 GeV; = 
-0.042 GeV; = 0.011/ GeV; Cl_ = -0.068 GeV; = -0.067,-0.095,-0.070 GeV, 
Afg = 1.90,-0.57,-0.11, and = -0.92,0.09,-0.21 for k = /, p and g respectively; = 
-0.062, -0.090, = -0.074, 0.094, and C^^ = 0.008, -0.02 for A; = g and p respectively. 
From Eq. ()17|) . we can extract isoscalar and isovector contributions, i.e. we consider the 

Tii- ^AA(in) ^isos (in) , ^isov (in) .,1 

splittmg L^, ' ' = ^ ' + hf.j^^ % with 

r-iisos (in) Ol^nB -^isov (in) jAj. f^k Oij^OinB (-\0\ 

- n^^XnB' ~ ^\ + nolnB + I^at' 

with = ak/{no + XnB)- 
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A choice of the isovector interaction in Eq. which is coherent with these results is 

where Gr is assumed to depend only on ms, whereas Gf and g are considered as functions 
of m4 only. The following vertices are deduced from it 



1-^202 ' / 1-^202 

In the MFA, they give rise to the following contributions to the nucleon self-energies 



S^isovAB X I^GrP sr^isovAB r ^3^f^ 

--Oar - Oab ^ . ja ' 



The isovector meson field equations can now be written as 

9S 



^ /^^isoijylA QY^isov AA \ ^ pj-^isov AA 

,0 



A 



dp 

In DDHFT j5|, |6|, |7| , the isoscalar vertices Tg and are defined by means of the relations 

Sf)=r,n,/m^, Sf) = r^n^/mL (19) 

where Ei*"\ E^t*"''' are the self-energies obtained with one boson exchange potentials in the 
Dirac-Brueckner-Hartree-Fock approach for symmetric nuclear matter. It must be noted 
that the right hand sides of these equations do not coincide in general with the dynamical 
self-energies. Therefore Eqs. fll9|) and p4|) - (fT3j) do not coincide, unless the rearrangement 
terms were omitted. In pjj, the author proposed that the vertices are solutions of someone 
of the differential equations ()14j] - ()l(ij] . instead. They must be solved together with the self- 
consistent condition for the meson fields. The explicit form of part of these solutions have 
been anticipated in so we summarize them below and add new results to the case (c) 
concerning isospin asymmetric nuclear matter. 

,rn \ 2 r^B 



a \ n 







nig J dn' 



f — ^ (Jr — ^ 

no + Xub 6 no + XriB p 




isov (in) 



j 

Jo 



dn' n' 



isov (in) 





isos (in) 



dT^vA 



dn'r 



isov (in) 



We have assumed that the input functions have been parameterized in terms of the partial 
nucleon densities ric, furthermore within the two last eqs. of case (c), the isospin parameter 
t is held constant in the integration. 

Thus we have obtained a set of relations, defining a hadronic field model suited to repro- 
duce the nucleon self-energies provided by other theoretical framework. In particular, the 
meson dependent vertices have been extended to deal with isospin asymmetric matter. The 
resulting vertices are given as functions of n^, and eventually S3, but they can be rewritten 
in terms of the meson fields or nucleon condensates. 

In the next section, these results will be compared with the standard DDHFT treatment, 
and the ability to adjust the nuclear matter phenomenology will be examined. 

III. RESULTS AND DISCUSSION 

A well known feature of isospin symmetric nuclear matter is its energy per particle, having 
a minimum at a baryonic density about uq = 0.16 /m~^. This property gives rise to bound 
states at zero temperature, and it is the least requirement that a model of nuclear matter 
should satisfy. Starting with Eq. Q, we have evaluated the energy- momentum tensor T^^" 
by the canonical procedure, and the energy per unit volume in the MFA is obtained by 
taking the in-medium expectation value of T°°. 



where we have introduced the Fermi momentum pp, related to the baryonic density by 
ub = 2p^p/{3tt^), and the effective nucleon mass M* = M — S^. The mean value of the 
isovector meson fields becomes zero for symmetric nuclear matter. 

Finally, the binding energy is defined as Eb = Emfa/^b — M, which should have a 




/ 1 

- + M*2 + F^n^cD + (Sr^ - F.a) + - ( 



mla'-mlu'), (20) 
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minimum value Eb — —16 MeV at the normal density no, to satisfy the nuclear matter 
phenomenology. Using the thermodynamical relation P = fiUB — Em fa we have obtained 
the pressure for the nuclear matter, with the chemical potential given by = Ep + 



being Ep = ^/pp + M*^. 

In first place, we have examined the effects of introducing the effective vertices by means 
of the algebraic Eqs. (fTIJ|) . or by solving a differential equation of the type listed in Eqs. 
flT ^ - (fT^ . It must be noted that the first option has the property of reproducing the energy 
density of the original Dirac-Brueckner calculations, whenever the parameterization of the 
vertices in terms of only the baryon number density is chosen. Otherwise, it is a definition 
without any special physical significance. The alternative procedure proposed in ^|, aims 
to make a connection between the hadronic vertices and another theoretical framework, by 
imposing the equality of the nucleon self-energies evaluated in both cases. 
For the purpose stated, we have chosen as input the results of . In this work, an ansatz for 
the vertices is proposed that fits several specific Dirac-Brueckner outcomes, but avoiding the 
unphysical behavior they exhibit in the zero density limit. The ansatze for the couplings 
and F^^ are rational functions of the relative baryonic density ns/nQ. The corresponding 
self-energies are obtained as si*"'* = T^^a, si*"'* = F^^cj, with the meson mean field values 
given by equations similar to those of Eq. (fT^ . For details, see Ref. jsj. 
In Fig. Qwe show the binding energy and pressure evaluated within the standard DDHFT 
treatment of [J]. Given the functions si*"'', si*"'' as above, we can also use them to deduce the 
vertices F^, Ty^ which reproduce these self-energies within the approach (a); i.e. solving Eqs. 
(I14|) . The corresponding equation of state is shown in the same figure. We have examined 
densities up to four times no, a range that may be reached inside neutron stars. It can be 
seen that similar results for the pressure are obtained for densities below 1.25no, but the 
binding energy shows appreciable differences. These results are qualitatively comparable, 
although for high densities the standard DDHFT treatment gives a softer growth for both 
Eb and P. These differences could be significative for those phenomena dominated by the 
regime of extreme densities, such as the structure of neutron stars. 

As the next subject we have examined the effect of multiple representations. Indeed, 
given a set of functions S^, depending on the macroscopic variables of nuclear matter, 
there exist an indefinite number of models like that of Eq. (j21), capable of reproduce them 
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within MFA. Some of these possibihties have been hsted as cases (a-c) in the previous sec- 
tion. The formal aspects have been summarized in Eqs. (fT^ - |TB|) . now we focus on some 
thermodynamical observables of symmetric nuclear matter. 

To avoid biased conclusions we have considered two inputs of diverse source. In addition to 
the results of |8j], based on one-boson exchange potentials, we have included the QCD sum 

g energy and the pressure, for each 
of these inputs under the three assumptions (a-c). We have found that despite the different 
parameterizations used for a given input, the results are practically undistinguishable in 
both cases. In Fig. El we show this feature for the calculations using the inputs obtained 

n 

from 1^], where it can be seen that the three curves are nearly coincident. The range of 
densities has been restricted to the region of validity of the QCD sum rules computations. 
The binding energy is monotonically decreasing, it seems to have a minimum near n ~ 2 no 
where the effective nucleon mass tends to zero. This description does not adjust to the 
nuclear matter phenomenology, however we keep Ref. jj] into consideration as it offers a 
physically meaningful input. It must be noted that the result (c) does not coincide with 
the preliminary result given in Ref. as the uj meson contribution to the energy was not 
properly taken into consideration there. 

A similar conclusion is obtained by using the self-energies extracted from O] , since the curves 
corresponding to the instances (b) and (c) (not shown in the figure) follow closely the solid 
curves of Fig. d in the full range < n < 4 no. 

The different behavior between the inputs obtained from 0| and j^, for the energy per 
particle, can be justified by examination of Fig. 01 In the upper panel the self-energies 
as functions of the density, show a monotonous increase in both cases. However the rate 
of growth of becomes lesser than the corresponding to for densities n > uq, in the 
parameterization of The QCD sum rules outputs instead, keeps a faster increase of the 
scalar as compared to the vector self-energy, over all the domain of densities. This behavior 
is emphasized in the bottom panel, where the difference — S^, is monotonous increasing 
in one case, but exhibits a extremum in the other one. Furthermore, a constant quotient 
is obtained for the QCD sum rules input, but it increases slowly for the other instance. 
Anyhow, the rational density dependence of Eq. (fTH|) for the iso-scalar self-energies could be 
retained in order to simulate the nuclear matter mechanism of saturation, for instance with 
as = —0.278, a^, = —2.93, and A = 0.697 we have obtained a binding energy of —14 MeV 
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Case 


as 


bs Cs 








Zs 




(a) 


0.805 


-0.170 0.350 


-0.082 0.386 


-0.032 0.148 


-0.024 


Si 


Si 


(b) 


1.264 


-0.232 0.294 


-0.031 0.413 


-0.018 0.096 


0.009 


S2 


Si 


(c) 


0.059 


0.000 0.066 


0.078 1.178 


-0.217 0.285 


-0.037 


1712 


nil 



TABLE I: Numerical values of the fitting coefficients for the density dependence of the vertices, 
obtained in the approaches (a-c) using the self-energies deduced from . 

at n ~ 1.1 riQ. These values must be contrasted with the parameterization deduced from jj], 
as = -2.309, = -0.66, and A = -0.261. 

The density dependence of the vertices, obtained by using the input of Ref. un- 
der the three approaches (a-c) is presented in Fig. HJ It must be noted that the ex- 
pressions for Ts^w become singular at zero density, nevertheless they have finite limits: 

— i> mldH^s"'^ / dn, — > 'm'^dT,v"'\o) / dn, for all the cases under consideration. 
A common feature of the vertices is a strong decrease up to n < 2 no, which moderates above 
this value. The cases (b) and (c) yield very similar results for both F^ and F^, differences 
with the case (a) may lead to qualitatively distinct descriptions for subsequent applications. 
All these results for the vertices Tk{z) can be summarized by an expression in terms of 
its natural variable z, = {z + a^)/ {hk + Ck{z + a) + dk{z + a)^), with k = s,w and the 
numerical coefficients are shown in Table 1. 

Another insight into the basic phenomenology of nuclear matter, can be given by the Lan- 
dau parameters. They can be evaluated as the Legendre projections of the second derivative 
of Emfa with respect to the baryonic density, for a deduction see jl^. These parameters 
are very useful in regard to collective phenomena of the dense nuclear environment, as for 
instance, phase transitions, the giant monopolar and quadrupolar modes, or the sound ve- 
locity [l^. In order to keep track of the momentum dependence, we assume a non-zero 
spatial component of the self-energy S^,, the baryonic current j and the omega meson tJ, 
taking the zero limit of these quantities at the end of the calculations. So, for example, we 
write Ep = \J {p — E)^ + M*^ into the integral of Eq. ()20|). Denoting by the occupation 
number of nucleon states with a well defined momentum p^, then we define 



kl 



dukdni 



(21) 
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M* dJ^s dT^v pk 9S„ 9E„ / dj Ug dT^y \ ^ duo duj 



Ei drii drii drii dn^ \ drii M* drii I ^ dn^ drii 

in the second equality, the hmit of isotropic nuclear matter has been taken. The derivative 
of the baryonic current can be evaluated in the MFA 

-H H = — dp 



drii E, drii' 6n^ Jo ^ (p2 + jvf* 2)3/2 ' 

The remaining derivatives are 

duj u dj d'Ey Pi S^, 
drii n drii ' drii Ei n + H'Ey' 

whose explicit form depends on the approach used (a-c). 

The Landau parameters fk are defined by means of the projections 

2m + l f\ p . 



27r 



1 



with p the cosine of the angle between pk and p;, and Pm the Legendre polynomial of order 
m. At the end of the calculations, the replacement \ pk \i \ Pi \ — > Pi? is made. By this 
procedure, we have obtained the Landau parameters of zero and first orders 

Ep dn dn ' 

It is costumary to normalize the Landau parameters with the density of states at the Fermi 
surface Sp = ^pfEf/t^"^, so that Fk = SFfk- 

For a given model of self-energies, the result for /q does not depend on the field parameter- 
ization of the vertices. The parameter fi instead, depends on uj which certainly differs for 
each of the instances (a-c). In Fig. El the density dependence of both Landau parameters is 
shown, significative departures for /i are found for densities above 0.25no. 
Consequently all physical magnitudes depending solely on /o, such as the nuclear compress- 
ibility K = 3pp{l + Fq)/Ef, do not distinguish among the field parameterization of the 
vertices. The same feature is shared by the giant monopole mode energy Em = chq/A^^^, 
with ao — 1-076 \/K/ jj,, and the first sound velocity vi = Pf\/ (1 + Fq)/{?,^Ef). On the 
other hand, the excitation energy of the quadrupole state is given by Eq = 02/ A^^^, with 
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^2 = PfV^/ (l-S/ii^l + Fi/3), therefore it is sensitive to the approach used. We have ob- 
tained in our calculations K ~ 440 MeV, ~ 145 MeV, and 02 — 77 MeV at the normal 
density, which must be compared with the empirical values 230 MeV < K < 270 MeV, 
uo = 80 MeV, and uj2 = 63 MeV. 

Another interesting phenomenological property that can be related to the Landau param- 
eters, are the zero sound modes. These are longitudinal collective modes propagating in 
nuclear matter at zero temperature, whose dispersion relation can be found as the zeros of 
the longitudinal dielectric function 

with Q = pqEp/{\ p I pf) and $ the Lindhard function. Col 
studied in detail within relativistic field models, see for instance 
lying collective modes can be classified as instability modes, in the low density regime, and 
zero sound modes with a characteristic linear dispersion relation pq (x\ p \. In Fig. |B[ 
the dispersion relation in terms of the baryonic density is shown for the standard DDHFT 
treatment and for the approaches (a-c), within the self-energy parametrization of [s]. These 
results are appropriate for the regime | p |— > 0, but keeping the linear dispersion relation. 
The instability and zero sound modes are clearly distinguishable, as the first one is closely 
related to the equation of state, all the corresponding curves practically coincide. The zero 
sound yields very different behaviors, the extremum values for its threshold are obtained 
within the approaches (b) and the standard DDHFT treatment. We have obtained n/uQ = 
1.25 and 2, respectively. In all the cases we have obtained two-folded zero sound, and taking 
into account that ^{Q) has non zero imaginary part for Q < 1, only the upper mode is 
undamped. 

IV. CONCLUSIONS 

We have examined in this work, two definitions for the vertices of a hadronic field model 
in terms of the self-energies provided by another theoretical framework. The first one is 
the standard DDHFT treatment, which uses Eqs. (fT^ to define the vertices F^, F^ be- 
tween nucleon and scalar and vector mesons, respectively. There, E^**^) denotes the nucleon 
self-energies, based on Dirac-Brueckner calculations with one-meson exchange potentials. 



ective modes have been 
Low 



3, 0,0,0, Q- 
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Although vector and scalar dependencies have been stated most practical applications 
have used only the first one, similar to our case (a). This particular choice has the property 
of reproducing the equation of state of the source calculations, but on the contrary the self- 
energies formally differ from the inputs S'^*"^ An alternative definition was proposed in 
it consisted in equating both, the input self-energy and that evaluated within the hadronic 
model with vertices depending on scalar combinations of hadronic fields. This procedure 
ensures a coherent overlap of the effective hadronic field model and a more fundamental 
theoretical description, at least in the baryon sector. Three different parametrizations of the 
vertices, cases (a-c), and two theoretical sources, namely the QCD sum rules of Ref. and 
the Dirac-Brueckner parametrization of Ref. , have been examined within this scheme. 
Thermo dynamical aspects of symmetric nuclear matter evaluated in MFA, are undistin- 
guishable for the instances (a-c) within the approach of see Eqs. ()14HTBj) . but they differ 
noticeably from the standard DDHFT results. Although the almost perfect coincidence 
within the (a-c) approaches for the energy and pressure, the vertices obtained present signi- 
ficative deviations in the medium to high densities regime. The results for Tg, are similar 
in the (b) and (c) instances, as they use parametrizations in terms of variables strongly re- 
lated in the MFA. We have given fittings for all the vertices in terms of its corresponding 
variables. 

In reference to the theoretical source for the inputs E*^*"), the parameterization given in 
yields observables in agreement with the nuclear matter phenomenology, as expected. The 
QCD sum rules instead, do not provide an acceptable description, at least in its present sta- 
tus. Taking the nucleon self-energies proposed in 4], we have presented formal expresions 
for the vertices of a hadronic field model for both isoscalar and isovector channels. However, 
examination of the binding energy obtained for symmetric nuclear matter, show a qualita- 
tive mismatching. This assertion is in opposition to the preliminary results of where 
the vector meson contribution to the energy density was not taken into account correctly. 
However we have shown that the functional dependence of the isoscalar self-energies given 

n 

in j4| is not completely reasonless, so as a redefinition of its numerical parameters produce 
a satisfactory equation of state. 

We have also considered in detail the Landau parameters for symmetric nuclear matter, 
under the four assumptions and the parametrization of ^. The hadronic model yields two 
non-zero Landau parameters Fq and Fi, the first one has a common behavior while the 
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second one strongly depends on the vertices definition. Therefore, nuclear observables de- 
pending on Fo only, such as isothermal compressibility, first sound velocity and the giant 
monopole energy, give almost coincident values. Those quantities depending on Fi may 

discriminate among the different approaches, for instance the zero sound modes exhibit very 
different threshold densities according to the choice made. 

To sum up, we have compared several formal schemes of merging external physical informa- 
tion into hadronic field models. For this purpose we have introduced interaction vertices 
and F^, which are functional of the hadronic fields, and we have determined them by else 
the standard DDHFT procedure or by requiring the reproduction of the self-energies used 
as input. Either of them, yield comparable results for the binding energy and pressure of 
symmetric nuclear matter. However, physical quantities can be found able to distinguish 
among the different approaches. 

QCD sum rules calculations for the nucleon self-energy, do not meet the requirements of the 
procedure, however further refinements could improve its performance. 
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100 




FIG. 1: The binding energy (top) and pressure (bottom) for symmetric nuclear matter at zero 
temperature, evaluated in the MFA using the standard DDHFT treatment and within the approach 
(a) . 
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G. 3: The nucleoli self-energies for symmetric nuclear matter, provided by references |a] and 
. In the upper panel the magnitudes of S^*"'* (solid lines) and H^"'^ (dashed lines), in the lower 
panel the difference (solid lines) and the quotient (dashed lines) between the scalar and vector 
components. 
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